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ABSTRACT 

We report new improved photometric redstarts of 1048 galaxies in the Hubble Deep Field (HDF). 
A standard % 2 minimizing method is applied to seven-color UB VIJH K photometry by Fernandez-Soto, 
Lanzetta, & Yahil (1999). We use 187 template SEDs representing a wide variety of morphology and age 
of observed galaxies based on a population synthesis model by Kodama & Arimoto (1997). We introduce 
two new recipes. First, the amount of the internal absorption is changed as a free parameter in the range 
of E(B — V) = 0.0 to 0.5 with an interval of 0.1. Second, the absorption due to intergalactic HI clouds is 
also changed by a factor of 0.5, 1.0, and 1.5 around the opacity given by Madau (1995). The total number 
of template SEDs is thus 187 x 6 x 3 = 3, 366, except for the redshift grid. The dispersion a z of our 
photometric redshifts with respect to spectroscopic redshifts is a z — 0.08 and 0.24 for z < 2 and z > 2, 
respectively, which are smaller than the corresponding values (a z = 0.09 and 0.45) by Fernandez-Soto 
et al. Improvement is significant, especially in z > 2. This is due to smaller systematic errors which are 
largely reduced mainly by including three opacities due to intergalactic HI clouds. We discuss redshift 
distribution N(z) and cosmic star formation rate based on our new photometric redshifts. 

Subject headings: galaxies: distances and redshifts — galaxies: photometry — galaxies: statistics 



1. INTRODUCTION 

Large redshift surveys of galaxies have played a quite 
important role for understanding the formation and evo- 
lution of galaxies. In 1990's, several new deep redshift sur- 
veys using large optical telescopes were carried out: e.g., 
Canada-France Redshift Survey (CFRS; Lilly et al. 1995), 
Hawaii Deep Field (Cowie et al. 1996), Autofib/LDSS 
survey (Ellis et al. 1996), the Canadian Network for Ob- 
servational Cosmology (CNOCl,2 surveys; Yee et al. 1996 
and Carlberg et al. 1998). Some of those surveys reach 
as deep as z ~ 1, and we have obtained some information 
about the formation and evolution of galaxies up to z ~ 1. 
However, little has been known about the properties of 
z > 1 galaxies from such redshift surveys due to the lim- 
ited depth of spectroscopic observations. Even if one uses 
LRIS with a Keck telescope, the limiting magnitude of 
spectroscopy is about Rn m ~ 25, while the limiting mag- 
nitude of photometry with Hubble Space Telescope (HST) 
WFPC2 is R Um ~ 28 - 29. 

There is, however, so called 'photometric redshift tech- 
nique', which enables us to measure redshifts of faint 
galaxies. Photometric redshift technique was originally 
suggested by Baum et al. (1963) to determine redshifts of 
some clusters of galaxies based on multi-band photometry 
using mean spectra of their member galaxies. There were, 
however, only a few applications of the technique (e.g., 
Koo 1985; Loh & Spillar 1986) before it revived recently to 
study distant galaxies, in particular, faint galaxies found in 
the Hubble Deep Field (HDF; Williams et al. 1996) (Con- 
nolly et al. 1995; Lanzetta, Yahil, & Fernandez-Soto 1996, 
1998; Gwyn & Hartwick 1996; Mobasher et al. 1996; Saw- 
icki, Lin, & Yee 1997; Connolly et al. 1997). Photometric 
redshift technique is the only method available at present 
to infer the redshifts of galaxies beyond the spectroscopic 
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limit. 

In this work we improve the accuracy of photometric 
redshifts of galaxies in the HDF by using a set of template 
spectral energy distribution (SED). We treat the inter- 
nal absorption of galaxies and the extinction due to in- 
tergalactic HI clouds more carefully than previous studies. 
We adopt a technique based on a standard x 2 minimizing 
method with a template set which consists of simulated 
SEDs. We apply it to a photometry catalog of galaxies in 
the HDF (Fernandez-Soto, Lanzetta, & Yahil 1999), which 
is the deepest catalog available with some spectroscopic 
redshifts. 

The detailed method, adopted parameters of simulated 
SEDs, and the comparison between our photometric red- 
shifts and spectroscopic redshifts are described in § 2. We 
show in § 3 the comparison of our improved photometric 
redshifts with previous studies and the validity of our re- 
sult. We present our improved photometric redshift for 
the HDF galaxies in § 4, where we also discuss the proper- 
ties of the HDF galaxies such as redshift distribution and 
cosmic star formation rate. Conclusions are given in § 5. 

2. METHOD 

We may classify photometric redshift techniques into 
two categories. One is the training-set method (Connolly 
et al. 1997, Wang et al. 1999) and the other is the \ 2 
minimizing method (Gwyn & Hartwick 1996; ; Mobasher 
et al. 1996, 1999; Sawicki et al. 1997; Giallongo et al. 
1998). 

The training-set method first makes an empirical for- 
mula by fitting polynominals to the photometry data of 
galaxies with known redshifts. It then applies the same 
formula to photometry data of new galaxy samples to de- 
termine their photometric redshifts. This method is supe- 
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rior to the other in the point that the more galaxies we 
observe for the training set, the more accurate photomet- 
ric redshifts we obtain. Therefore, Yee (1998) called it 
the empirical training-set method. For low-z galaxies, the 
training-set method seems to work more robust than the 
X 2 method, because it is based on real observed galaxy 
spectra and because the empirical formula can be deter- 
mined accurately using a large number of spectroscopic 
samples (Connolly et al. 1997). Recently, comparably 
accurate redshifts are obtained by a simpler training-set 
method which adopts binning of galaxies according to 
observed colors (Wang, Bahcall & Turner 1998; Wang, 
Turner, & Bahcall 1999). Csabai, Connolly, & Szalay 
(1999) also succeeded in obtaining accurate photometric 
redshifts by combining the training-set method and the \ 2 
method. However, it is difficult to make a good training 
set for high rcdshift galaxies, since the number of galaxies 
with spectroscopic redshifts is still small. Also we should 
be cautious in applying the empirical formula obtained 
at low redshifts to high redshift galaxies, since significant 
evolution is expected in their SEDs. 

Hence, in this study, we use the x 2 minimizing method 
which can include spectral evolution of galaxies directly. 
This method searches for the best-fit SED and the redshift 
of a galaxy by comparing the observed SED of the galaxy 
with templates prepared in advance. In general, templates 
are made either from observed spectra of nearby galaxies 
or from simulated spectra based on a population synthe- 
sis model. Coleman, Wu, & Weedman (1980; hereafter 
CWW) and Kinney et al. (1993, 1996) are widely used 
sources of observed spectra for the HDF galaxies. Lanzetta 
et al. (1996, 1998) and Fernandez-Soto et al. (1999) used 
CWW's and Kinney et al.'s spectra. However, the use of 
the observed SEDs of nearby galaxies may make it dif- 
ficult to directly take the effects of galaxy evolution into 
account. Accordingly, in this paper we use simulated spec- 
tra based on a recent population synthesis model. Gwyn 
& Hartwick (1996) also used simulated spectra based on 
Bruzual & Chariot (1993). Sawicki et al. (1997) made 
their templates based on both observed spectra (CWW) 
and simulated spectra by Bruzual & Chariot (1996; here- 
after BC96). 

2.1. Template SEDs 

We use a stellar population synthesis model by Ko- 
dama & Arimoto (1997; hereafter KA97) to make tem- 
plate SEDs. KA97 includes the stellar evolutionary tracks 
after the asymptotic giant branch in addition to the con- 
ventional tracks. Hence the model predicts UV flux of 
galaxies reasonably well, or at least, as good as any previ- 
ous models. KA97 was successfully used to obtain photo- 
metric redshifts of low redshift galaxies (Kodama, Bell, & 
Bower 1999). Our work is the first to use KA97 to obtain 
photometric redshifts of high-z galaxies. 

The template SEDs consist of the spectra of pure disks, 
pure bulges, and composites made by interpolating the 
two as shown in Table |l|. The parameters for the SEDs 
are the power-law index a; IMF of the initial mass function 
(IMF), the time scale of star formation tsf, the time scale 
of gas infall from a galactic halo into a disk Ti n f a n, and 
the time when the galactic wind blows £gw- The star for- 
mation rate of a galaxy is set to be zero after £gw- For 
disks, we adopt x lMF — 1.35, tsf = 5Gyr, Ti n f a n = 5Gyr, 



and £qw = 20Gyr (i.e., longer than the present age of the 
universe), which are close to the values estimated for the 
disk of our Galaxy. We adopt from KA97 a; IMF = 1.10, 
r SF = O.lGyr, T infa n = O.lGyr, and t GW = 0.353Gyr for 
bulges, which are known to reproduce the average color of 
ellipticals in clusters of galaxies. We make intermediate 
SED types by combining a disk component and a bulge 
component with the same age. The ratio of the bulge lu- 
minosity to the total luminosity in the B band, which we 
define as B/T, is changed from 0.1 to 0.9 with an interval 
of 0.1. Pure disk SEDs correspond to young or active 
star-forming galaxies, and pure bulge SEDs correspond 
to elliptical galaxies. We also prepare very blue SEDs of 
age<lGyr, corresponding to blue star-forming galaxies re- 
ported by recent deep surveys. In total, our basic template 
set consists of 187 SEDs. 

2.2. The Absorption Effects 

To simulate the observed SEDs, we take into account 
two absorption effects. One is the internal absorption due 
to dust in each galaxy. More light is scattered and ab- 
sorbed by dust grains in shorter wavelengths. Therefore, 
internal absorption is one of the key parameters for photo- 
metric redshift technique which is sensitive to wavelength 
dependent effects. For simplicity, we assume that all galax- 
ies have the same extinction curve, although we change the 
absolute amount of extinction, indicated by E(B — V), as 
a free parameter. 

Typical features of extinction curves of nearby galaxies 
including our Galaxy are as follows: (1) The amount of ab- 
sorption basically increases from the infra-red to the ultra- 
violet. (2) Some extinction curves have a bump at 2175A, 
which is usually explained mainly by graphite. The exis- 
tence of the bump is reported for our Galaxy (Seaton 1979; 
Cardelli et al. 1989), LMC (Fitzpatrick 1985), and M31 
(Bianchi et al. 1996). However, no bump is reported in 
the extinction curve of SMC (Prevot et al. 1984; Bouchct 
et al. 1985). Gordon, Calzetti, & Witt (1997) showed that 
whether the bump is present or not cannot be explained 
by scattering or geometrical effects, and that galaxies with 
high star formation activity tend to have no bump in their 
extinction curves. 

There are some studies in which an analytical formula 
was fitted to observed extinction curves (Cardelli et al. 
1989; Fitzpatrick 1986,1998). Calzetti (1997a) proposed 
an analytical formula for nearby star-forming galaxies, in- 
cluding the effect of scattering through star-forming re- 
gions. He inferred that his formula holds for distant young 
galaxies because these galaxies are probably similar to 
nearby star-forming galaxies. 

We calculated photometric redshifts using the Milky- 
Way-like extinction curve by Cardelli, Clayton, & Mathis 
(1989) and the SMC-like curve by Calzetti (1997b). As a 
result, Calzetti's curve gave better results than Cardelli et 
al.'s. Therefore, we have decided to adopt the curve by 
Calzetti (1997b), 

F obs w = Fvmo-° AE{B - v)kW , (1) 

where F b s (X) is the observed flux and Fq(X) is the intrin- 
sic flux without extinction. If we normalize the absorption 
strength by k(B) — k(V) = 1, fc(A) is expressed as: 

k(X) = 2.656(-2.156+ 1.509/A - 0.198/A 2 + 0.011/A 3 ) 
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+ 4.88 
= [(1.86 



for 0.12 fxm < A < 0.63 /zm 
' 1.73 



-0.48/A)/A-0.1]/A- 

for 0.63 fim < A < 1.0 /jm (2) 



It is not clear whether or not most of distant galaxies have 
a SMC-like extinction curve because of the ambiguity of 
the UV flux. 

Figure [l] shows the observed extinction curves and the 
model curve calculated by Calzetti (1997b). The curves 
are normalized so that E(B — V) is equal to 1. We make 
the template SEDs using the Calzetti's extinction curve 
with E(B -V) = 0.0 (dustless), 0.1, 0.2, 0.3, 0.4, and 0.5. 
Figure g shows the ratio of the flux after extinction with 
respect to the intrinsic flux as a function of wavelength for 
three values of E(B — V). 

We also include the effect of absorption due to the inter- 
galactic HI clouds. On the assumption that the fr-value, 
which is the indicator of the velocity dispersion of the 
clouds, is 35km/s, Madau (1995) computed the median 
optical depth of the line blanketing as 
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2.3. x 2 -Fitting and Constraints 

The fitting procedure to obtain reduced \ 2 is summa- 
rized as follows. We prepare a set of original 187 SEDs 
based on KA97, which cover a wide range of the spectral- 
type of observed galaxies. Then, we incorporate the ab- 
sorption effects due to the internal dust and the intergalac- 
tic HI clouds (§ 2.2). Next, we make redshifted SEDs from 
z — 0.05 to z = 8.0 with an interval of 0.05. Finally, each 
SED at a given redshift is convolved by the band responses 
used for the observation. We compute the flux of template 
galaxies observed in the i-th band by 
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where fl e ™ plate d\ is the flux received in [A, A + g?A] and R\ 
is the band response function of the z-th band. 

We compare the convolved flux Tj with the observed flux 
Fi in each band and calculate \ 2 by: 
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where j = 1,2,3, and 4 correspond to Ly-a, Ly-/3, Ly-7, 
and Ly-5 absorptions, respectively, A.,- = (1216, 1026, 973, 
950A), and Aj = (3.6 x 10~ 3 , 1.7 x 10" 3 , 1.2 x 10" 3 , 9.3 x 
10~ 4 ). The amount of absorption of lines at shorter wave- 
lengths is negligible compared with these Ly-a to Ly-5 
lines and Lyman continuum. Madau (1995) also showed 
that the continuum absorption below 912A including Ly- 
man limit systems is expressed with a 5% accuracy by: 
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where x c = 1 



' Zc = iff and Xem = 1 + z em- Here, A obs is 
the wavelength of the observer's frame and Al=912A and 
z em is a redshift of an emitting source. The total opacity 
of intergalactic HI clouds, T e ff, is expressed by the sum of 
the opacity due to line absorption r l J^ e and that due to 
continuum absorption • : 
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(5) 



The relative flux due to intergalactic HI absorption is 
shown in Figure |^. We also estimate the amount of the 
HI absorption integrated within standard broad-band sys- 
tems as shown in Figure f|. 

We include a variation in opacity as follows. In addition 
to Madau's original median opacity ir e g ) , we also use the 
opacity of 0.5T e ff(= f~g) and l.hr e ff{= t^). This is be- 
cause we should take into account a variation of Lyman 
absorptions statistically. Here, and r~g are within 
±1(7 of the opacity derived by Madau (1995). As we will 
see later, using three different values around the median 
value turns out to be more preferable to using just the 
median opacity. 



vation, Cj is an observational error in the i-th band. If 
flux in the i-th band is not detected, Ft is set to be zero. 
a denotes the normalization factor which minimizes the 
reduced \ 2 ■ For a given set of Tj, the value of a which 
minimizes \ 2 



2 is calculated by: 
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Then, we find the best-fit SED which gives 'the minimum 
reduced % 2 ' out of all the templates, and the redshift of 
the best-fit SED is adopted as the photometric redshift of 
the galaxy. 

The most important point to get reliable photometric 
redshifts is the validity of the template SEDs. However, 
even with an ideal set of templates, large estimation errors 
in Zphoto caused by, for example, small signal-to-noise ra- 
tios (S/N) of photometry are inevitable in some cases. We 
should carefully check the results obtained by our simple 
fitting method. In this work, we impose two constraints 
in order to get robust and accurate estimations. First, we 
reject the best-fit SEDs that give too bright or too faint 
absolute magnitudes. We then remove from our sample 
the best- fit SEDs which have 'minimum reduced x 2 ' larger 
than some threshold. It is found that the photometry data 
of too low S/N tend to give a large value of 'minimum re- 
duced x 2 '. 

3. VERIFICATION OF THE METHOD 

3.1. The Spectroscopic Sample 

We apply our method to galaxies in the HDF, where 
the deepest images of field galaxies were taken. We use 
the photometric catalog by Fernandez-Soto et al. (1999), 
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which is based on seven-band photometry, i.e., UBVI with 
the HST (Williams et al. 1996) and JHK taken with the 
4m Mayall telescope at the Kitt Peak National Obser- 
vatory (Dickinson et al. 1999). This catalog also gives 
one of the most recent estimates of photometric redshifts 
(hereafter, z p hoto) of the HDF galaxies. This catalog in- 
cludes total magnitude and photometric errors for galaxies 
of Iab < 28 with a detection threshold of /i/ AB =26.1 [mag 
arcsec" 2 ], though the magnitude limit is Iab — 26 in the 
outskirts of the field. We summarize the limiting magni- 
tudes in Table ||. The number of galaxies in the catalog 
is 1067 in total, and 946 galaxies are in the good S/N re- 
gion, while 121 galaxies are in the outskirts. Out of them, 
108 galaxies have spectroscopic redshifts (hereafter, z spe c) 
(Cohen et al. 1996; Steidel et al. 1996; Lowenthal et al. 
1997; Zepf, Moustakas, & Davis 1997; Cohen et al. 2000), 
which we hereafter call the spectroscopic sample. 

3.2. Photometric Redshifts and SED Parameters 

We first measure z p hoto for the spectroscopic sample, us- 
ing UB VIJHK seven-band magnitudes by the standard \ 2 
fitting of the template SEDs. Figure |5| shows the correla- 
tion diagram between our z p hoto and z spGC . We cannot find 
the best-fit SED and fail to measure z p hoto for two galax- 
ies, and one galaxy at z ~ 3 has a very large estimation 
error, which is called the catastrophic error by Sawicki et 
al. (1997). As seen in Fig||, the redshifts of four galaxies 
at z spec ~ 2 — 3 are underestimated by Az ~ 0.5. A similar 
trend was seen in Fernandez-Soto et al. (1999) and some 
other previous studies (Gwyn & Hartwick 1996; Lanzetta 
et al. 1996, 1998; Sawicki et al. 1997). However, we can 
identify such problem atic galaxies by the additional con- 
straints explained in § 2.3 and reject them from the sample 



to reduce systematic errors. 

We examine the minimum reduced \ 2 value, which is 
an statistical indicator of the mean residual over all the 
observed bands. Figure ^| shows the reduced \ 2 values 
against the differences between z p hoto and z spcc . Figure^ 
demonstrates that it is difficult to reject the galaxies with 
bad z p hoto based on the minimum reduced % 2 . However, 
we find that the x 2 value in the B band, Xbj is a good 
indicator for removing bad z p hoto estimates. We exclude 
17 galaxies with \b > 40 from the sample. The break- 
down number of these galaxies for different redshift bins 
is shown in Table [| 

It is probably inappropriate to pursue 100% success in 
z p hoto, especially for a sample of very distant galaxies. It 
is inevitable that some fraction of the sample fails to give 
correct z p hoto- The important point is how to identify such 
false z p hoto and reject them from the sample so that the 
final sample has the least contamination even with less 
completeness. 

We now have photometric redshifts for 91 out of 108 
galaxies in the spectroscopic sample, which are shown in 
Figure 0. In Figure 0, the galaxy pointed by an arrow has 
very poor S/N. There are two galaxies with large devia- 
tions at z spcc ~ 0.5. Hogg et al. (1998) noted that the 
spectroscopic redshifts of these galaxies are determined by 
just one emission line, suggesting that the line identifica- 
tion is possibly wrong. Table ^presents the accuracy of our 
z p hoto estimates for the 91 galaxies. We present in Table || 
the standard deviations for three cases, total 91 galaxies, 
90 galaxies without the z spcc = 2.93 galaxy having a catas- 



trophic error, and 86 galaxies with |z spec — z p hoto| _ 0.5. It 
is found that our method gives almost the same accuracy 
as previous estimates for galaxies at z < 2, and gives more 
accurate z p hoto for galaxies at z > 2. 

In order to show that the \ 2 minimizing method selects 
a reasonable SED from the SED templates for each of the 
observed galaxies, we present in Figure || '2-D scatter di- 
agrams' of parameters for the best-fit SEDs. The area of 
the filled circle is proportional to the number of SEDs on 
the grid. In the panels, r = 1 corresponds to the Madau 
value of T e ff and r = — and + represent r~g = 0.5r eJ y 

and T~^g = 1.5r eJ y, respectively. In panel (c), all galaxies 
at z < 1 are fitted with r = 1, but this is an artifact. 
For galaxies at z < 1, the intergalactic Lyman absorption 
does not enter any of the seven passbands, and r = 1 is 
assigned to such galaxies. 

Figure || demonstrates that the distributions of param- 
eters seem to be broadly consistent with the properties of 
actual galaxies. For example, panel (a), the diagram of 
T3/T vs. E(B — V)', shows that the population with large 
E(B — V) is dominated by disk-dominated galaxies. Also 
in panel (g), the 'z spec vs. age' diagram, no old galaxies 
are seen at high redshifts. We conclude that Figure || gives 
indirect support to the validity of our method. 

Finally, we examine the effects of different S/N on the 
accuracy of z p hoto using simulations. We prepare 1,077 
input galaxies by a random selection from total 541,926 
grids of our SED template, which consist of 3,366 SEDs 
and 161 redshifts. These input galaxies cover a wide range 
of spectral types. We add to them Gaussian noise so that 
photometric errors in UB VIJHK bandpasses be A = 10%, 
which is roughly the same as those of the HDF galaxies. 
We obtain z p hoto of these input galaxies in exactly the 
same way as for the HDF galaxies. The results are shown 
in Figure || (the case of A = 20% is also shown for com- 
parison) and the values of a z are summarized in Table |[ 
The best-fit SED for a input galaxy is not always the same 
as input SED itself because we add Gaussian noise to the 
magnitudes of input galaxies. However, if we do not add 
noise, exactly the same SED as input is answered. 

4. COMPARISON WITH PREVIOUS WORK 

Hogg et al. (1998) reported a comparison among z p hoto 
by various authors for the same low-z (z < 1.4) galaxies in 
the HDF. We summarize their results (Table 3 of Hogg et 
al. 1998) again in Table ||, where we also give our result 
for exactly the same sample for reference. It is found that 
the accuracy of our z p hoto estimate for these low-z galax- 
ies is comparable to those by the other authors. However, 
it should be noted that the authors given in Hogg et al. 

(1998) estimated their z p hoto as a completely blind test, 
while our estimate is not blind. We have optimized some 
of the parameters in our method on the basis of the pho- 
tometry and spectroscopy catalog by Fernandez-Soto et al. 

(1999) which already includes most of the galaxies used for 
Hogg et al.'s blind test. In that point, our result can not 
be directly compared with those by the other authors. The 
optimization in our method is, however, done using z > 2 
galaxies only so that the accuracy of z p hoto for these high- 
z galaxies be better. Accordingly, the accuracy of z p hoto 
for such low-z galaxies is not affected significantly by the 
optimization. 
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We compare our z p hoto of high-z galaxies in the HDF 
with those by Fernandez-Soto et al. (1999), which are 
based on the same sample and the same photometry as we 
use. Therefore, this is a direct comparison of the meth- 
ods. Fernandez-Soto et al. had two galaxies with catas- 
trophic errors, while wc have one. The galaxy for which 
both of the authors give a catastrophic error is located 
at z spcc = 2.93. This galaxy is fitted by a young dusty 
spectrum of 0.5Gyr and E(B — V) = 0.5, and cannot be 
rejected by our x 2 threshold. The other galaxy which has 
a catastrophic error in Fernandez-Soto's study is rejected 
by the % 2 threshold in our work, though it is also fitted by 
a dusty young disk SED. Table [7] summarizes the z p hoto 
accuracy for galaxies of z > 2, with catastrophic errors 
excluded. Our work gives a much improved accuracy for 
high- z galaxies. Wang, Turner, & Bahcall (1999) achieved 
er z = 0.08 and 0.30 for z < 2 and 2 < z < 4 galaxies, 
respectively, which are comparable to our results. They 
have adopted a training-set method which uses colors. 

Figure [10| compares z p hoto among authors (data of Gwyn 
& Hartwick (1 996) and Sawicki et al. (1997) arc take n 
from Sawicki [ http : / / www. astro .utoronto . ca / "sawicki / 1 ) . 
The results by Fernandez-Soto et al. (1999) and ours have 
been newly added. The results by Sawicki et al. and 
Gwyn & Hartwick are based on the optical 4-band pho- 
tometry for 74 galaxies with spectroscopic redshifts, while 
Fernandez-Soto et al.'s and our work are based on the 7- 
band photometry for 108 galaxies. 

Fernandez-Soto ct al. and we obtain more accurate 
Zphoto over all the redshift range than Gwyn et al. and 
Sawicki et al., presumably because 7-band photometry has 
longer baseline in wavelength and hence is more robust for 
estimating z p hoto than 4-band photometry. The optical 4- 
band photometry is worse at z = 1 — 2 because in that red- 
shift range the 4000A break of galactic spectra drops out 
of the optical bands and no significant feature is available. 
The photometric redshifts by Gwyn & Hartwick seem to 
be worse than Sawicki et al.'s. This is probably because 
Sawicki et al. included Lyman absorption into their tem- 
plate SEDs while Gwyn et al. did not. 

Now we examine the absorption effects on z p hoto esti- 
mation. In Gwyn et al.'s estimation, there are two galax- 
ies with catastrophic errors at z spGC — 3. Sawicki et al. 
(1997) suggested that these catastrophic errors can be im- 
proved by including the reddening effects into the template 
SEDs and thereby reducing ' aliasing', the misidentifica- 
tion between a 4000A break and other breaks in shorter 
wavelengths. We try our method using the template SEDs 
without internal absorption. The results are shown in Fig- 
ure [y]. Figure [ll] indicates that the absence of internal 
absorption in template SEDs increases catastrophic errors 
at Zgpcc — 3. This is probably because 'no dust SEDs' 
have steeper 4000 A breaks than real observed galaxies, 
and therefore it is difficult to distinguish Lyman breaks 
(or some other breaks) from 4000A breaks with no dust 
SEDs. We have adopted dusty SEDs in this study. The 
values of E(B — V) are taken in the range of 0.0 (no dust) 
to 0.5. These values are consistent with those for Lyman- 
break galaxies at z = 2 — 3, E(B — V) ~ 0.3, obtained 
by Sawicki & Yee (1998). Calzetti & Heckman (1999) also 
suggested E(B — V) < 0.5 for galaxies at z > 2 — 3 by 
calculating the cosmic star formation rate and E{B — V) 



values of galaxies iteratively. 

In Figure |io|(a)-(c), we see that galaxies at z spcc = 2 — 4, 
especially for the seven-band photometry, tend to be gen- 
erally estimated at lower redshifts than z spcc by the three 
previous studies. We are able to remove this systematic 
error by including three values of opacity due to the in- 
tergalactic Lyman absorption around the median opacity 
derived from Madau (1995): We include T e ff, the exact 
value by Madau, T eff{= 0-5r ej (j), and t^(= 1.5r eJ y). The 
systematic errors at z — 2 — 4 emerge if we remove tJ^ and 
r~g and use only T e g (Figure [l2|). It is clear that includ- 
ing a variation of the opacity is essential. We also examine 
the template SEDs with 3r e ff, 4r eJ y, 0.3r e ff, and 0.25r ej fj, 
and find that the best combination is (r e ff i T eff) ■ We 
should point out that the ambiguities of the stellar syn- 
thesis model used for the template SEDs, especially those 
in the UV flux, may be coupled with the ambiguity of the 
internal absorption and the statistical fluctuation of the 
intergalactic Lyman absorption. If we look at Figure ||(c) , 
best-fit t values seem to change with redshift: is more 
favorable for z spcc = 2 — 4 and r~g is more favorable for 
Zspec> 4. However, we do not insist that this trend is real 
considering the fact that we do not have much knowledge 
about the UV flux of galaxies. Nonetheless, we can obtain 
more accurate 2 p hoto at 2 < z < 4 by carefully treating 
the intergalactic Lyman absorption. 

In summary, we find that the internal absorption is ef- 
fective in reducing the catastrophic errors, while the inter- 
galactic absorption is effective in reducing the systematic 
error at z — 2 — 4. 

5. PROPERTIES OF THE HDF GALAXIES BASED ON THE 
IMPROVED PHOTOMETRIC REDSHIFTS 

5.1. How Many Template SEDs? 

We have used 187 template SEDs covering a wide range 
of spectral types. Generally speaking, accurate z p hoto, 
even for high-z galaxies, can be obtained using either ob- 
served spectra (CWW; Kinney et al. 1993) or simulated 
spectra (KA97; BC96), though absorption effects (internal 
and intergalactic absorption) must be carefully treated. If 
the purpose of the analysis is just to obtain z p hoto (e.g., 
Fernandez-Soto et al. 1999), a few template SEDs would 
be sufficient. However, if we want to determine accurately 
the SED shape of a target galaxy in order to derive useful 
information on the SED, much more template SEDs are 
necessary as we will see below. 

Figure |l^ shows z p hoto of the HDF galaxies using four 
SEDs from CWW as template SEDs, which is almost the 
same procedure as Fernandez-Soto et al. adopted (Note 
that we reproduce the same trend in Figure l|| as seen 
in Figure |l0|(c)). The accuracy of z p hoto based on the 
four SEDs is worse than that of our best estimate (Fig- 
ure 0(c)), but the differ ence is within a factor of about 
2. However, there is a large difference in the minimum 
reduced \ 2 ■ O ur minimum reduced x 2 values are plot- 
ted in Figure [| Figure 14 is the same plot as Figure || 
but for the four-SED case. The minimum reduced % 2 val- 
ues obtained with the four-SED templates are much larger 
than those with 3,366 (187 x 6E(B - V) x 3r) templates. 
The distribution of reduced % 2 values can be a measure 
to determine how many template SEDs are necessary and 
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sufficient for a given purpose. 

5.2. Properties of the HDF Galaxies 

We obtain photometric redshifts of 925 g alax ies out of 
the 946 galaxies in the good S/N region (See § |3.l|). We call 
these 925 galaxies the photometric sample. The redshift 
distribution N(z) of the photometric sample is shown in 
Figure [15[ For comparison, the results by Fernandez-Soto 
et al. are superposed. 

Sawicki et al. (1997) argued that N{z) could be an 
indicator of the error in photometric redshift technique. 
They claimed that aliasing among spectral breaks (4000A, 
Balmer, 2800A, and 2635A) easily leads to an unrealistic 
N(z), which has typically two prominent peaks at z ~ — 1 
and z ~ 2. The peaks are likely caused by a template set 
which does not include the internal or the intergalactic 
absorption. Our N(z) does not have such two prominent 
peaks. Instead, our N(z) shows a single peak at z ~ 0.5 — 1 
and a moderate decrease at z > 1 at all magnitude bins. 

Fernandez-Soto et al. found that galaxies at z = 
2 — 4 have systematically lower photometric redshfits (Fig- 
ure |0|(c)). We find in this paper that this is attributed to a 
systematic error which is sensitive to the intergalactic ab- 
sorption. This 'small aliasing' due to the systematic error 
would make an unlialistic peak at z ~ 2 in N(z), especially 
for faint galaxies. For 26 < I < 28, the N(z) of Fernandez- 
Soto et al. has a prominent peak at z ~ 2. This is probably 
due to the small aliasing of galaxies at z > 2, because the 
number of Fernandez-Soto et al.'s galaxies with z > 2 is 
smaller than that of our galaxies. 

Next, we examine E(B — V) of galaxies. Figure |l6| shows 
the histograms of E(B — V) of the best-fit SEDs in the five 
redshift ranges. The arrow in each redshift range indicates 
the mean value. We find that the E(B — V) value of galax- 
ies remains almost constant (E(B — V) = 0.13 — 0.18) 
for z = — 6, though it has a small peak at z = 3 — 4 
bin. These values are roughly consistent with those of 
Lyman-break galaxies at z — 2.5 — 3.5 derived by Saw- 
icki & Yee (1998; E(B - V) ~ 0.3), because our pho- 
tometric sample probably includes not only star-forming 
Lyman-break galaxies but also passively evolving ellipti- 
cal galaxies. Calzetti & Heckman (1999) iteratively calcu- 
lated the cosmic star formation rate (SFR) and the mean 
E(B — V) of galaxies using the Galactic conversion fac- 
tor between the gas column density N(H) and E(B — V), 
i.e., E(B - V) = N(H)/5.9 x 10 21 (Bohlin, Savage, & 
Drake 1978). The vertical lines plotted in Figure lq are 
predictions of E(B — V) based on their two typical mod- 



els. The model 'A' corresponds to the cosmic SFR which 
remains almost constant from z — 1.5 to 5. The model 'B' 
corresponds to the cosmic SFR which reaches its peak at 
z = 1 — 2 and monotonically decreases at z = 2 — 5. Our 
results are found to be more consistent with the model 'A'. 

6. SUMMARY AND CONCLUSIONS 

We improve photometric redshifts of 1048 galaxies in the 
HDF. A standard x 2 minimizing method is adopted for ob- 
taining Zphoto- We use 187 template SEDs representing a 
wide variety of morphology and age of observed galaxies 
based on a population synthesis model by Kodama & Ari- 
moto (1997) with two new recipes. First, the amount of 
the internal absorption is changed as a free parameter in 
the range of E(B — V) = 0.0 to 0.5 with an interval of 
0.1. Second, the absorption due to intergalactic HI clouds 
is also changed by a factor of 0.5, 1.0, and 1.5 around the 
opacity given by Madau (1995). The total number of tem- 
plate SEDs is thus 187 x 6 x 3 = 3, 366, except for the 
redshift grid. The introduction of the internal absorption 
is found to be effective in reducing the catastrophic errors 
in Zphoto while the absorption due to the intergalactic HI 
clouds helps to improve the systematic error in z p hoto at 
z = 2 — 4. The dispersion o~ z of our photometric redshifts 
with respect to spectroscopic redshifts is a z = 0.08 and 
0.24 for z < 2 and z > 2, respectively, which are smaller 
than the corresponding values (a z = 0.09 and 0.45) by 
Fernandez-Soto et al. (1999). Significant improvements 
are obtained, especially at high-z. A comparison with pre- 
vious work and the properties of the best-fit SEDs for the 
spectroscopic sample verify the validity of our photometric 
redshifts. 

The redshift distribution of all the 925 galaxies in the 
photometric sample shows a peak at z ~ 1 and a moderate 
decrease at z > 1 without any other peak. The E[B — V) 
value of galaxies remains almost constant for z = — 6 
with a possible weak peak at z — 3 — 4, which is consistent 
with the model of a constant cosmic SFR given by Calzetti 
& Heckman (1999). 

We would like to express our gratitude to T. Kodama & 
N. Arimoto for kindly allowing us to use their new popula- 
tion synthesis code and to anonymous referee for construc- 
tive comments. H.F. wishes to thank the Japan Society 
for the Promotion of Science for a financial support. This 
work is supported in part by Grants-in-Aid (07CE2002, 
11640228, 10440062) from the Ministry of Education, Sci- 
ence, Sports and Culture of Japan. 
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Table 1 

Parameters of template SEDs adopted in this work 



SED type 


x IMF(a) 


T SF [Gyr] 


T infaii [Gyr] 


taw [Gyr]( b ) 


Age [Gyr] 


pure disk 
(37kinds) 


1.35 


5.0 


5.0 


20.0 


0.010, 0.013, 0.016, 0.020, 0.025, 
0.032, 0.040, 0.050, 0.063, 0.790, 
0.100, 0.126, 0.158, 0.200, 0.251, 
0.316, 0.398, 0.501, 0.631, 0.794, 
1.0, 1.259, 1.585, 2.0, 3.0, 4.0, 
5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 
12.0, 13.0, 14.0, 15.0 


pure bulge 
(15kinds) 


1.10 


0.1 


0.1 


0.353 


I. 0, 2.0, 3.0, 4.0, 5.0, 
6.0, 7.0, 8.0, 9.0, 10.0, 

II. 0, 12.0, 13.0, 14.0, 15.0 


Bulge+Disk^ 
(9xl5kinds) 










I. 0, 2.0, 3.0, 4.0, 5.0, 
6.0, 7.0, 8.0, 9.0, 10.0, 

II. 0, 12.0, 13.0, 14.0, 15.0 



( a ) Power-law index of initial mass function. 

'k-^GW = 20 means that the galactic wind does not blow until the present epoch. 
( c )B/T=0. 1-0.9, where B and T are the bulge and the total luminosity. 



Table 2 
The photometry of the HDF 



Band 


Exposure time 


PSF(FWHM) 


limiting magnitude 




[hours] 


[arcsec] 


[mag] 


F300W( U) 


42.7 




26.98(AB) for 10<r limit 


F450W(B) 


33.5 




27.86(AB) for 10<r limit 


F606W( V) 


30.3 


0.12 


28.21(AB) for 10<r limit 


F814W(/) 


34.3 




27.60(AB) for 10(7 limit 


J 


11.0 


1.0 


23.45(STD) for 5cr limit 


H 


11.3 


1.0 


22.29(STD) for 5cr limit 


K 


22.9 


1.0 


21.92(STD) for 5a limit 



Note.— Upper 4 rows: HST-WFPC2 (UBVI) (Wiiliams ct al. 1996) and 
lower 3 rows: KPNO-IRIM (JHK) (Dickinson et al. 1999). 



Table 3 

The breakdown number of galaxies in the spectroscopic sample 



•^spcc 


0.0 - 6.0 


0.0 - 2.0 


2.0 - 6.0 


Total with z S pec 


108 


79 


29 


Good[ x | < 40] 
Bad[ X | > 40] 


91(84.3) 


66(83.5) 


25(86.2) 


17(15.7) 


13(16.5) 


4(13.8) 



Note. — The values in the brackets mean the ratio(%) 
to the total number within each rcdshift bin. 
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Table 4 

Summary of accuracy of our z pho to for the spectroscopic 

SAMPLE 



z spcc Number of galaxies u z 







1 


91 


0.329 


0.0 - 


6.0 


2 


90 


0.171 






3 


86 


0.101 






1 


66 


0.139 


0.0 - 


2.0 


2 66 


0.139 








64 


0.081 






1 


25 


0.585 


2.0 - 


6.0 


2 


24 


0.239 






3 


22 


0.145 



[1] For all galaxies. 

[2] One galaxy with a catastrophic error is removed. 
[3] Galaxies of z S pcc- Zphotol > 0.5 are removed. 



Table 5 

Summary of accuracy of our z P hoto for the simulation with 10% Gaussian noise 



Zspcc 


0.0 - 6.0 


0.0 


- 2.0 


2.0 


- 6.0 


Number of input galaxies 


1077 (1071) 


272 


(269) 


805 


(802) 




0.102 (0.081) 


0.157 


(0.116) 


0.076 


(0.066) 



Note. — Numbers in the parenthesis are for those with |z spC c- z p hoto < 0.5. 



Table 6 

Comparison of accuracies for galaxies with <z spcc < 1.4 according to Hogg et al. (1998) 



Author 


Number of galaxies 


Band 


Fraction 1 ' 






Az < 0.1 [%] 


Az < 0.3[%] 


Gwyn 


23 


UBVI 


65.2 


91.3 


Mobashcr 


20 


UBVI 


35.0 


80.0 


Sawicki-4band 


20 


UBVI 


65.0 


85.0 


Sawicki-7band 


20 


UBVIJHK 


70.0 


95.0 


Connolly 


20 


UBVIJ 


65.0 


95.0 


Fcrnandcz-Soto 


19 


UBVIJHK 


73.7 


94.7 


This Work 2 ) 


19 


UBVIJHK 


78.9 


94.7 



^Fraction (%) of galaxies which have Az < 0.1 and 0.3 is given for different authors. 

2 )This work is not a blind test as explained in the test. However, the same sample is 
used as in Fernandcz-Soto's result. 
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Table 7 

Comparison of a z of high-z (z > 2) galaxies between our estimates and Fernandez-Soto et al.'s 







a z (Fernandez-Soto) 


a z (This work) 


Sample 


A (N=27) 


0.58 


0.29 


Sample 


B (N=26) 


0.45 


0.30 


Sample 


C (N=24) 


0.42 


0.24 



Note. — Numbers in the parenthesis are number of galaxies in 
each sample. [A]: A galaxy which has a catastrophic error in our 
estimate is rejected. [B]: Two galaxies which have a catastrophic 
error in Fernandez-Soto et al's estimate are rejected. [CI: Five 
galaxies which have either a catastrophic error or a large \ value 
in our estimate are rejected. 
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Fig. 1. — The extinction curves of typical galaxies. The vertical axis means the absorption strength A(A), which is normalized by 
E(B — V) = 1. The curve given in Calzctti (1997b) is indicated as the solid line labeled 'Model'. It is noted that Milky Way and LMC have 
the 2175A bump while SMC and Calzetti's curves do not. 
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Fig. 2. — The relative flux after absorption by interstellar dust in a galaxy. We assume Calzetti's extinction curve. 
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Fig. 4. — The absorption strength in magnitude due to intergalactic HI clouds in case of observing a galaxy at z = 3.5 with broad band 
filters. The U, B, and V bands here means F300W, F450W, and F606W niters of the HST, respectively. 
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Fig. 6. — Minimum reduced \ 2 plotted against the error in rcdshift estimation for all galaxies in the spectroscopic sample. The horizontal 
axis indicates the differences between z p hoto and z spC c in the logarithmic scale. 
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Fig. 8. — The parameter distributions of the best-fit SEDs for the spectroscopic sample. The area of the filled circles is proportional to the 
number of SEDs in the grid. For example, in panel (d) (B/T vs. r), the grid of (B/T,t)=(0.0, 1) includes 41 galaxies, the largest circle of 
all. See the text for details. 
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Fig. 9. — Output .Zphoto versus input 2 spe c of the simulation. The left and right panels are for the 10% and 20% input noises in flux, 
respectively. 
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Fig. 10. — Comparison of 2 p hoto versus z Bp c C among various authors. Results by Sawicki ct al. and Gwyn et al. are based on optical 4-band 
photometry for 74 galaxies with spectroscopic redshifts. Results by Fcrnandcz-Soto et al. and this work are based on 7-band photometry for 
108 galaxies. 
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Fig. 12. — z p hoto based on the template SEDs using just one value, T e g, for the opacity of intergalactic Lyman absorption. 
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Fig. 13. — %,hoto obtained based on CWW's four SEDs as a template set. 
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Fig. 14. — Minimum reduced \ 2 obtained with the CWW's four SEDs as a template set. Compare with Figure ^ 
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N(z) of the photometric sample. The solid and dotted lines indicates our results and Fernandcz-Soto ct al.'s. 
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Fig. 16. — Histograms of E(B — V) for five redshift bins. Arrows indicate the mean of E(B — V) in each redshift bin. Vertical lines in each 
bin denote the predictions by Calzetti & Heckman (1999). The dashed lines are for model 'A' and the dotted lines are for model 'B' (See text 
for details). 



